Islands in the Landscape 



T. CliftonH, Andrei Linde@ and Navin Sivanandam§ 

Department of Physics, Stanford University, Stanford, CA 94305 



Abstract 

The string theory landscape consists of many metastable de Sitter vacua, populated 
by eternal inflation. Tunneling between these vacua gives rise to a dynamical system, 
which asymptotically settles down to an equilibrium state. We investigate the effects of 
sinks to anti-de Sitter space, and show how their existence can change probabilities in 
the landscape. Sinks can disturb the thermal occupation numbers that would otherwise 
exist in the landscape and may cause regions that were previously in thermal contact to 
be divided into separate, thermally isolated islands. 
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1 Introduction 



One of the most intriguing features of string theory is its prediction of a multitude of vacuum 
states [H [2l El IH [5] Stabilizing these states [6] and coupling the resulting embarrassment of 
riches with a population mechanism [8] provided by eternal inflation [9], [10] gives rise to the 
string theory landscape [llj . Physics in the landscape can be both rich and perplexing. We 
aim to elucidate some relevant ideas for understanding this physics, with the hope of improving 
our comprehension of the multiverse. 

Inflationary expansion divides the universe into many exponentially large domains, each 
corresponding to different metastable vacuum states. In this picture tunneling between dif- 
ferent vacuum states causes bubbles of new vacuum to be continually nucleated. Those with 
positive vacuum energy are initially static, but soon accelerate in their expansion until the 
velocity of their walls asymptotically approaches that of light. If the conditions are right inside 
these bubbles then a stage of slow-roll inflation will occur and the resulting observers will see 
themselves in an infinitely extended, open Friedmann universe. Percolation of successive bub- 
bles inside of each other give us a universe that is eternally inflating and constantly producing 
new inflationary universes, where structure can form, and life can evolve. 

The first step in a complete understanding of this scenario must be to find out which vacua 
are possible in string theory, and to describe their typical properties |3]. Once these vacua 
have been identified we will then need to study cosmological evolution during eternal inflation 
in order to determine the global structure of the universe [8]. During inflation the number of 
horizon-sized dS regions of space-time is continually, and exponentially, increasing. These dS 
regions then 'populate' the many possible vacua of string theory, realizing the great variety of 
the theory in a diverse and eternal universe. 

The resulting picture is incredibly complex. Ultimately, our goal is to explain the properties 
of our part of the multiverse, and to predict the results of future observations. To achieve 
this goal one needs to calculate the probabilities of various outcomes in an eternally inflating 
multiverse. This is a thorny problem, and a subject of much contention. Studying the global 
structure of an eternally inflating spacetime leads to comparisons of infinite volumes, and 
hence a consequent dependence on cutoff procedures [8j [121 H31 EEU EES! EES! EH]- Escaping 
cutoff problems is not impossible if one considers individual observers and concentrates on 
their individual histories, ignoring the rest of the universe [T%1 \W\ |2"D"1 |2"TI |2"2"] ; here, however, 
we must face the problems of initial conditions and are, perhaps, led to worry about Euclidean 
quantum gravity and the wave-function of the universe [231 1211 [25]. More importantly, this 
description tends to miss some of the important features of eternal inflation. 

In this paper we will leave anthropic considerations aside and concentrate on other prop- 
erties of the string theory landscape. We will focus on the existence, or otherwise, of thermal 
equilibrium between populations of dS vacua. As we will see, under certain conditions, a system 
of dS vacua described in comoving coordinates settles down to a state in which the popula- 
tions of these vacua are in thermal equilibrium with one another. More precisely, the ratios 
of comoving volume occupied by one vacuum or another will depend on the exponential of 
the entropy difference between them. An important limitation of this simple picture is that 
it is valid only in comoving coordinates, which do not reward different rates of cosmological 
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expansion in different parts of the universe. Nevertheless, the picture of many dS universes in 
a state of thermal equilibrium is very simple and intuitively appealing, and therefore it can be 
very useful for understanding various features of the string theory landscape. 

On the other hand, this simple picture may be invalid when the landscape has sinks (ter- 
minal vacua which can be tunneled to, but not from). In particular, in [26J [27] it was shown 
that for a simple system consisting of 2 dS vacua and one AdS sink, naive expectations of 
thermal equilibrium are incorrect if the decay rate to the sink is sufficiently fast. Since one 
expects sinks to be common in the landscape it may be the case that the disruption of 
thermal equilibrium between metastable dS vacua is a generic feature, and the string theory 
landscape may consist of many thermally isolated 'islands.' The goal of this work is to elucidate 
this possibility and to investigate in more detail the situations in which the usual thermally 
equilibrium populations are disturbed. We will find explicit solutions for a variety of simple 
configurations that may occur in the landscape, and use these results to form a picture of how 
the vacua of a more realistic landscape may be populated. 

It will be found that the presence of sinks in the landscape can significantly alter the 
dynamics of the inflating multiverse. One of the most dramatic and unexpected examples of 
this is that when a number of high energy vacua decay to a single lower energy vacuum, which 
can decay to a sink, the probability fluxes soon become dominated by the slowest decaying, 
most stable vacuum. In the limiting case of this vacuum being completely stable it makes no 
contribution at all to probability fluxes, as these are the results of tunneling events between 
metastable vacua (see below). However, if the smallest chance of tunneling out of this vacuum 
is allowed, we unexpectedly find that this tiny current comes to be the dominant source of the 
probability flux. This slowest decaying vacuum may then remain out of thermal contact with 
other vacua, whilst all faster decaying vacua eventually approach thermal equilibrium with each 
other. This 'tortoise and the hare' scenario shows explicitly the non-trivial effect of sinks on 
the dynamics of inflation in the string theory landscape: They may lead to the existence of 
thermally isolated, slowly decaying vacua while all other, more rapidly decaying vacua are left 
in thermal equilibrium. 

We will also investigate the possibility that some of the AdS or Minkowski sinks could 
potentially act as impassable barriers between systems of dS vacua, thus carving the landscape 
into totally disconnected 'islands'. Such a situation would result in different regions of the 
multiverse being completely isolated from one another, whilst maintaining thermal equilibrium 
internally. We argue that the large number of vacua and dimensions in the landscape, coupled 
with the 'vacuum dynamics' we find to exist in the presence of sinks, makes the existence of 
such isolated regions improbable, though perhaps not impossible. 

In section [2] we will review the basic mechanisms of tunneling between vacua. Following this, 
in section [3] we will outline the results of [251 EI], f° r a simple landscape of two dS vacua and one 
AdS sink. Section 0] discusses the possibility of sinks dividing the landscape into disconnected 
islands. In section Owe shall summarize the results of our investigations into more complicated 
toy landscapes, highlighting some of the counter- intuitive features that emerge. In section [6] we 
summarize our results. Mathematical details can be found in the appendix. 
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2 Tunneling in the Landscape 



There are two related mechanisms for making transitions between vacua: one due to tunneling 
[28J and another due to stochastic diffusion processes [HI [29]. A somewhat more detailed 
discussion of these mechanisms, and the issues associated with them, can be found in [27] . We 
summarize the salient points below. 

Tunneling between vacua produces bubbles of new vacuum, that look like infinite open 
Friedmann universes to observers inside. If the tunneling goes to dS space, then the bubble 
expands exponentially with the velocity of its walls approaching that of light. (In comoving 
coordinates these bubbles approach some maximal value and freeze. This maximal value de- 
pends on the time when the bubble is formed, and is exponentially smaller for bubbles formed 
later on [30].) If the tunneling goes to a state with a negative vacuum energy V, the infinite 
universe inside it collapses within a time of the order |V|~ 1//2 , in Planck units. These negative 
energy, AdS vacua then play the role of sinks for probability currents in the landscape. 




Figure 1: Coleman- De Luccia tunneling may go in both directions. 

Let us consider two dS vacua, dSj, with vacuum energy density V{ = V(4>i), Fig. [TJ Without 
taking gravity into account, the tunneling may go only from the upper minimum to the lower 
minimum, but in the presence of gravity tunneling may occur in both directions, which is 
emphasized in Fig. [TJ According to Coleman and De Luccia [28], the tunneling probability 
from dSi to dS2 is given by 

r 12 = e- B = e~ s ^ +s \ (1) 

where S((j)) is the Euclidean action for the tunneling trajectory, and Si = S(0i) is the Euclidean 
action for the initial configuration = 0i, 

* = (2) 
This action has a simple sign-reversal relation to the entropy of de Sitter space, Si: 

Si = -Si = +^ • (3) 
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Therefore the decay time of the metastable dS vacuum tdecay ~ T 12 1 can be represented in the 
following way: 

+ , - e sW+Si _ t S&) (A) 

''decay ° b r c • V / 

Here t r ~ e Sl is the so-called recurrence time for the vacuum dSi. 

Whereas the theory of tunneling developed in [28] was quite general, all examples of tun- 
neling studied there described the thin-wall approximation, where the tunneling occurs from 
one minimum of the potential and proceeds directly to another minimum. In the cases where 
the thin-wall approximation is not valid, the tunneling occurs not from the minimum but from 
the wall, which makes interpretation of this process in terms of the decay of the initial vacuum 
less trivial. 

The situation becomes especially confusing when the potential is very flat on the way from 
one minimum to another, V" < V, in Planck units. In this case the Coleman-De Luccia (CDL) 
instantons describing decay of a dS space do not exist [21]; they become replaced by Hawking- 
Moss (HM) instantons. According to Hawking and Moss [31], the probability of tunneling from 
the minimum 1 to the minimum 2 is then given by 



12 = e -Sto P +Si = exp / - " + " ) . (5) 



24tt 2 24?r 

The HM instanton is described by the Euclidean version of dS space corresponding to the top 
of the potential barrier, (ft = <p top . 

Unlike the thin-wall CDL solution, the HM solution does not interpolate between the two 
different minima of V(4>), and therefore debates on the validity of the HM result continue even 
now [32]. One may wonder why we should consider such instantons instead of considering the 
instantons corresponding to the dS space in the next minimum; the resulting tunneling action 
would be much smaller. Moreover, one may consider a string theory landscape with many 
minima and maxima separated by a sequence of barriers. Then one could wonder whether the 
HM tunneling suppression applies only to the tunneling between the nearby vacua, or if it can 
describe direct tunneling to distant minima, ignoring all intermediate barrier except the last 
one [3, [32]. One of the best attempts to clarify this situation was made by Gen and Sasaki [33], 
who described the tunneling using Hamiltonian methods in quantum cosmology, which avoided 
many ambiguities of the Euclidean approach. But even their investigation does not allow us to 
answer the last of these questions. 

A proper interpretation of the Hawking-Moss tunneling was achieved only after the de- 
velopment of the stochastic approach to inflation [SJ, Q21 |2H1 [31] . One may consider quantum 
fluctuations of a light scalar field <ft with m? = V" <C H 2 = V/3. During each time interval 
St = H~ l this scalar field experiences quantum jumps with the wavelength ~ H~ l , and with a 
typical amplitude S<p = H/2tt. As a result, quantum fluctuations lead to a local change in ampli- 
tude of the field 0, which looks homogeneous on the horizon scale if -1 . From the point of view 
of a local observer, this process looks like a Brownian motion of the homogeneous scalar field. 
If the potential has a dS minimum at <\>\ ^> H 2 /m, then eventually the probability distribution 
to find the field with the value at a given point becomes (almost) time-independent, 

( 24tt 2 24vr 2 \ 
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The distribution P(4>) gives the probability to find the field at a given point, and has 
a simple interpretation as the fraction of comoving volume of the universe in each of the dS 
vacua, or, equivalently, a fraction of time the field spends in a vicinity of its value <fi along 
the Brownian trajectory. Up to a sub-exponential factor, this distribution shows the density 
of points with a given value of the field <fi along its Brownian trajectory. This implies that, 
up to a sub-exponential factor, the typical time required for the field, at any given point in 
comoving coordinates, to move from its equilibrium value <pi and climb to the top of the barrier is 

proportional to P _1 (0 top ) ~ exp ^ — \ [7J. Once the scalar field climbs to the top of 

the barrier, it can fall from it to the next minimum, which completes the process of "tunneling" 
in this regime. That is why the probability to gradually climb to the local maximum of the 
potential at <fi = top and then fall to another dS minimum is given by the Hawking-Moss 
expression [HI EH [7J, EE]. It is also why tunneling to distant minima separated by many 
barriers is accomplished by a sequence of transitions from one minimum to another nearby 
minimum, rather than by one big jump. This last statement does not follow from the Hawking- 
Moss derivation of their result, but is apparent from the stochastic approach to inflation. 

A necessary condition for the derivation of Eq. ([6]) using the stochastic approach to inflation 
in [H [HI [29j EI] is the requirement that m? = V" <C H 2 = V/3. This requirement is satisfied 
during slow-roll inflation, but it is violated for all known scalar fields at the present (post- 
inflationary) stage of the evolution of the universe. Thus the situation with the interpretation 
of the Coleman-De Luccia tunneling for V" > V/3 is somewhat unsatisfactory. However, since 
the validity of the Coleman-De Luccia approach was confirmed at least in some limiting cases 
(in the absence of gravity, and in the slow- roll regime discussed above), in this paper we will 
follow the standard lore, assume that this approach is correct, and study its consequences. 

Following [35] (see also [TTJ ESI EH]), we will look for the probability distribution P to find 
a given point in a state with vacuum energy V,, and will try to generalize the results for the 
probability distribution obtained above by the stochastic approach to inflation. The main idea 
is to consider CDL tunneling between two dS vacua, with vacuum energies V\ and Vz, such that 
V\ < V2, and to study the possibility of tunneling in both directions, from V\ to V 2 , or vice 
versa. 

The action on the tunneling trajectory, S ((/>), does not depend on the direction in which the 
tunneling occurs, but the tunneling probability does depend on it. It is given by e~ s ^ +Sl on 
the way up, and by e ~ s ^ +S2 on the way down [35] (for a recent discussion of related subjects 
see also [37]). Let us assume that the universe is in a stationary state, such that the comoving 
volume of the parts of the universe going upwards is balanced by the comoving volume of the 
parts going down. This can be expressed by the detailed balance equation 



p e SW+Si _ p 2 e -S{4>)+S 2 



(7) 



which yields (compare with Eq. 



©) 





independently of the tunneling action S(4>). 
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This probability distribution also has a nice thermodynamic interpretation in terms of dS 
entropy S [39] : 

§ = e*-si = e AS . (9) 
Pi 

Here, as before, Si = —S{. This result and its thermodynamic interpretation have played a 
substantial role in the discussion of the string theory landscape 

Following [261 EZ] we shall now show this notion of thermal equilibrium is a fragile one, and 
can be destroyed by a simple extension of our model to include sinks. 



3 A Toy Landscape: 2 dS and 1 AdS 

Stationarity of the probability distribution (Q was achieved because the lowest dS state did 
not have anywhere further to fall. Meanwhile, in string theory all dS states are metastable, so 
it is always possible for a dS vacuum to decay [5]. Further, it is important that if it decays by 
the production of bubbles of 10D Minkowski space, or by production of bubbles containing a 
collapsing open universe with a negative cosmo logical constant, then the standard mechanism 
of returning back to the original dS state no longer operates. Therefore Minkowski vacua, as 
well as AdS vacua, work like sinks for the flow of probability in the landscape. Because of the 
existence of these sinks (also known as terminal vacua), the fraction of the comoving volume 
in the dS vacua will decrease in time. 

Although decays from one SUSY AdS vacua to another are forbidden, uplifting [6] breaks 
super symmetry. Uplifted dS vacua can then decay to AdS by the formation of bubbles of 
collapsing universes. According to [26] , the typical decay rate for this process can be estimated 

as T ~ exp (CMp/m\j^\ . For a gravitino mass, m 3 / 2 , in the 1 TeV range one finds suppression 

in the range of V ~ 10~ 1034 [26], which is much greater than the expected rate of decay to 
Minkowski vacua, or to higher dS vacua, which (in vacua like ours) is typically suppressed by 
factors of order 10~ 10 . Other possible decay channels for the uplifted dS space were discussed 
in [101I1I1I12]. 

Let us consider a simple model describing two dS minima and one AdS minimum (denoted 
by 1, 2, and S in Fig. [2]). Here, as in the rest of this paper, we work in comoving co-ordinates. 
To get a visual understanding of the process of bubble formation in comoving coordinates, one 
may paint black all of the parts corresponding to one of the two dS states, and paint white 
the parts in the other dS state. Then, in the absence of sinks in the landscape, the multiverse 
will become populated by white and black bubbles of all possible sizes. Asymptotically, it will 
approach a stationary regime - on average becoming gray, with the level of gray becoming 
asymptotically constant. Suppose now that some parts of the universe may tunnel to a state 
with a negative cosmo logical constant. These parts will collapse, so they will not return to the 
initial dS vacua. If we paint such parts red, then the universe, instead of reaching a constant 
shade of gray, eventually will look completely red. 

To describe this process, instead of the detailed balance equation (j7j) one should use the 
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Figure 2: A potential with two dS minima and a sink. 

"vacuum dynamics" equations [201 EE] : 

A = -J ls - J 12 + J 21 , (10) 

Pi = -J2s - hi + 3X2 • (11) 

Here Jij = Pj T^, where is the decay rate of the vacuum j to bubbles of vacuum i. In 
particular, J ls = Pi e~ Cl is the probability current from the lower dS vacuum to the sink, i.e. 
to a collapsing universe, or to a Minkowski vacuum, J 2s = P2 e~° 2 is the probability current 
from the upper dS vacuum to the sink, Ji 2 = P\ e~ Sl+ l 5 W is the probability current from the 
lower dS vacuum to the upper dS vacuum, and J21 = P2 e~ S2+ ' 5( ^ is the probability current 
from the upper dS vacuum to the lower dS vacuum. Combining this all together, gives us the 
following set of equations for the probability distributions: 

Pi = -Pi (Vis + r 12 ) + P 2 T 2 i , (12) 
P 2 = -P 2 (r 2s + r 21 ) + p 1 r 12 . (13) 

We ignore here possible sub-exponential corrections, which appear, e.g., due to the difference 
in the initial size of the bubbles etc. 

Because of the decay to the sink, Pi and P 2 gradually become exponentially small. But this 
does not mean that the whole universe goes to the sink: The physical volume of white and black 
parts of the universe continues growing exponentially, in the regime of eternal inflation. One of 
the ways to account for this growth is to slice the universe by hypersurfaces of time t measured 
in units of H~ x . In this case, volume of all parts of the universe during time At = 1 grows e 3 
times. One can describe this effect by adding the terms 3p to the r.h.s. of Eqs. ( |T2l . ( JTBl [27J. 
As a result, the functions Pj describing the total volume of the different dS vacua will grow 
exponentially even in the presence of the sink (if the rate of decay to the sink is not too large). 
In this case, the functions P, will correspond to the 'pseudo-comoving' probability distribution 
[27J. In this paper we will be interested only in the ratios of the volumes of dS spaces, Pi/Pj- 
These ratios are not affected by the overall growth of all parts of the universe, and therefore 
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the ratios Pi/Pj are the same for the comoving and pseudo-comoving probability distributions. 
Therefore for simplicity we will not add the terms 3Pj to the r.h.s. of our equations, i.e. we 
will study comoving probabilities. 

To analyze the solutions of equations (fT2|) and (fl3|) . let us first understand the relations 
between their parameters. Since entropy of dS space is inversely proportional to energy density, 
the entropy of the lower level is highest, Si > S2. As the tunneling is exponentially suppressed, 
we have S 2 > \S((f>)\, so we obtain a hierarchy Si > S 2 > |S(0)|, and therefore T 12 <C r 21 <C 1. 
We will often associate the lower vacuum with our present vacuum state, where Si ~ 10 120 . 

For simplicity, we will study here the possibility that only the lower vacuum can tunnel to 
the sink, T 2s = 0. The two equations (IT2|) and f[T5j) can then be solved exactly to give the 
general solution 

~ A 



p 1 = c ie -!(ri2+r ls +r 2l) (t-*o) cogh 



fit - k 



(14) 



where C\ and t are constants and A 2 = (r 12 + T ls + T 21 ) 2 — 4r ls r 21 . Substituting back into 
the Pi equation then gives the ratio 

P2 r 12 - r 21 + Tu + Atanh [j(t - t )] 

Pi ~ 2T 21 
Tu - r 21 + T ls + A 

~* ^ (15) 

as t — > 00. These solutions show us that in a comoving coordinate system, in the presence of a 
sink, both of Pi and P 2 are decaying, whilst their ratio approaches a constant value (i.e. their 
rates of decay become equal). 

One may consider two interesting regimes, providing two very different types of solution. 
Suppose first that Ti s <C T 21 , i.e. the probability to fall to the sink from the lower vacuum is 
smaller than the probability of the decay of the upper vacuum. In this case one recovers the 
result obtained without the sink in the previous section: 

£ = ^ = e S2 - Sl « I- (16) 
r\ 1 21 

It is interesting that this thermal equilibrium is maintained even in the presence of a sink if 
Tis <C r 2 i- Note that the required condition for thermal equilibrium is not Ti s Ti 2 , as one 
could naively expect, but rather Ti s <^ T 2 i. We will call such sinks narrow. 

Now let us consider the opposite regime, and assume that the decay rate of the uplifted dS 
vacuum to the sink is relatively large, Ti s 3> T 2 i, which automatically means that Ti s ^> r i2 . 
In this "wide sink" regime the solution of Eq. ffl5|) is 

^ _ Eii _ gSa-ISWI-d w e S a -|5(0)| > j nj\ 

Pi r 2 i 

i.e. one has an inverted probability distribution. This result has a simple interpretation: if the 
"thermal exchange" between the two dS vacua occurs very slowly as compared to the rate of 
the decay of the lower dS vacuum, then the main fraction of the volume of the dS vacua will 
be in the state with the higher energy density, because everything that flows to the lower level 
rapidly falls to the sink. 
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4 Is the Landscape Transversable? 



Before discussing the details of how AdS sinks can effect thermal structure in the landscape, 
we will make a note of another scenario that can arise in their presence: The separation of the 
landscape into disconnected regions, or islands. 

Consider the simple potential shown in Fig. [3j One may naively expect that all energetically 
favorable transitions between vacua should occur; however, this may not be the case in the 
presence of sinks. In the absence of gravity we know that quantum mechanical tunneling 
can only occur from higher energy vacua to lower energy ones. In such a picture all vacua 
slowly tunnel to lower and lower energy, until they eventually end up in the lowest energy 
ground state. Including the effects of gravity changes this situation by allowing the possibility 
of tunneling upwards, to higher energy. Although going upwards is less probable than going 
down, between any two vacua, the system soon settles down into a situation where the fluxes 
going up and down are equal. This behavior occurs when the population in the lower energy 
state far exceeds the number in the high energy state, compensating the unlikeliness of jumping 
upwards by increasing the number of vacua that could potentially make this transition. This 
process is known as recycling, and, for a system of dS vacua only, allows the whole system to 
approach thermal equilibrium, where the ratios of occupation numbers of any two vacua are 
given by the equilibrium ratios flUD- 

AdS and Minkowski sinks, which can be tunneled to, but not back from, spoil this recycling 
mechanism and can disrupt the global thermal structure that would otherwise exist. One 
manifestation of this effect is the altered occupation fractions Pi/Pj, discussed in the previous 
section; another is the possibility of these sinks separating the landscape into isolated regions, 
which may be in thermal equilibrium internally, but not with each other. We call such regions 
'islands'. 




S 



Figure 3: A potential with two dS minima seperated by a sink. 

Fig. [3] shows a simple situation which may occur in the landscape; two dS vacua with 
an AdS sink between them. If tunneling is allowed between the two dS vacua, then thermal 



10 



equilibrium between them can occur. Conversely, if tunneling is not allowed then these two 
vacua will be totally isolated from one another. 

Consider first the possibility of a Coleman-De Luccia transition from vacuum 1 to vacuum 
2. As was shown in [31], the CDL instantons exist only for V" > V, where tunneling is usually 
assumed to occur from one side of a potential barrier to the other, as in Fig. [TJ The scalar 
field then rolls down the potential to the minimum, where reheating may occur if conditions are 
right. For the situation shown in Fig. [31 however, such rolling will be down to an AdS minimum, 
which corresponds to a collapsing open universe. One may speculate about the possibility of a 
tunneling from a collapsing space, or about the re-emergence of the universe after the collapse, 
but we do not know how to study this regime in a controllable way@ Alternatively, one may try 
to find the CDL instantons describing a direct tunneling from vacuum 1 to vacuum 2, jumping 
over the AdS minimum. 

On the other hand, if we consider V" <C V, then the CDL instantons do not exist. In this 
case we still have the Hawking-Moss instanton, and its interpretation in stochastic inflation. 
However, according to the stochastic interpretation of the HM tunneling, there can be no 
tunneling directly from 1 to 2. Recall that we view this scenario as representing the scalar field 
tunneling to the top of the barrier through a sequence of quantum jumps and subsequently 
rolling down to an adjacent vacuum. Clearly this is not a plausible transition between 1 and 2, 
since the only adjacent vacuum is the central AdS sink. 

Thus, in general one may encounter situations where it is not possible to travel from one 
part of the landscape to another. However, it may happen that whereas the probability of the 
transitions discussed above may be extremely small, they cannot be strictly forbidden; our 'no 
trespassing' result may be just a consequence of the approximation used to study such processes. 
Moreover, considering a mult i- dimensional potential with many moduli and fluxes will also add 
complexity, as tunneling between vacua 1 and 2 could occur in the other dimensions. With an 
increasing number of dimensions it becomes increasingly likely that tunneling will be able to 
occur either directly between 1 and 2, or via some intermediary vacuum. Despite these caveats, 
we find it interesting that the existence of sinks allows for the possibility of the landscape being 
separated into isolated islands. 

In the following section we will consider the possibility of tunneling events directly between 
two islands, separated by a sink. Such events could be due to a direct tunneling via the 
Coleman-De Luccia instantons in the regime V" > V, or due to a sequence of the tunneling 
events though intermediate vacua. We will show that even if these events are allowed and the 
islands are not totally isolated, they may still be thermally isolated from one another as the 
sink can disrupt the equilibrium that would otherwise exist. However, our analysis leads us to 
believe that if the number of dimensions in the landscape is large enough then the existence of 
sinks will not be able to disrupt the tunneling in all of these dimensions, and thermal contact 
will therefore be maintained. Furthermore, the 'tortoise and the hare' behavior we find indicates 
that if a single high-energy vacuum is allowed to decay to vacua on each of the islands then this 
vacuum will act as a bridge, restoring thermal contact. For a landscape with very many vacua 
and dimensions, it therefore seems implausible that thermally isolated regions should occur. 
We discuss this in more detail below. 

4 We are grateful to Tom Banks for a discussion of this issue. 
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5 Non- Trivial Thermal Structure in the Landscape 



In this section we will present our results and analysis of more general situations than the simple 
one shown in Fig. [2j We explore the situations in which thermal equilibrium (which is to say, 
detailed balance) between vacua is disturbed. We hope that this discussion will be of interest 
to a wider audience, and can be understood without having to follow lengthy calculations. The 
interested reader can then proceed to the appendix, where the mathematical minutiae are given. 

5.1 The General Equation 

For a system of n vacua and a single sink we can describe the evolution of the probability 
measures, P, by: 

dP 

= ~ E T ^ + E T a p i - r * p <- ( 18 ) 

This equation is completely general: Single sinks to which many vacua can decay with differ- 
ing amplitudes, and multiple sinks to which individual vacua can decay (again with differing 
amplitudes) amount to the same thing. The exact form of the general solution to this set of 
equations is found in Appendix A.l: 

Pi = ^ ] C%j e j > 
j 

where the cy are constants of integration, and the rrij are constants formed from the transition 
rates, T. All Pj have the same functional form, and their coefficients are related by factors 
which are functions of the Ts only. 

At late times the ratios between different vacua then asymptote to the constant values: 

Pi 

rij = lim — , (19) 

t— >oo Pj 

which can easily be found, using the solutions for Pj. Of course, finding solutions in a generic 
landscapes is not easy, so we will restrict ourselves to a number of simple examples. 

In section [3] we found the solution for the potential illustrated in Fig. [2J and proceeded to 
find the asymptotic limit of the ratio P%jP\ as t — > oo. We now find simple estimates for the 
asymptotic form of the ratios Pi/Pj in more generic situations. 

5.2 An Extended ID Potential 

Consider now the extended one-dimensional potential in Fig. HJ Here we allow transitions 
between neighboring de Sitter minima in both directions, and transitions from 1 to the AdS sink. 
The equations determining the vacuum dynamics, and their solutions, are given in Appendix 
A.2.1. 

For definiteness, we will assume that V 3 > V 2 > V%, and therefore r 32 > T 2 3, and T 2 i > r 12 . 
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Figure 4: A one- dimensional potential with three positive and one negative minima 
Firstly, when Ti s <C r 2 i, we find: 

P\ r 2 i p 3 r 23 

— ~ — » 1 and — ~ — < 1. 

r 2 J- 12 r 2 1 32 

In this limit the thermodynamic ratios © are maintained because T 2 i T 12 and T 32 ^> r 23 . 
When T ls ^> T 21 and T 21 <C r 32 , we find: 

Pi 

— ~ — < 1 and — ~ — < 1. 

r 2 1 Is "1 1 32 

Here the thermal ratio of Pi to P 2 is broken, whilst that of P 2 to P 3 is maintained. Lastly, 
when Ti s ^> r 2i and r 2i ^> T 32 , we find: 

P\ r 2 i P 3 r 2 i 

— ~ — < 1 and — ~ — > 1. 

-T2 1 Is 2 1 32 

Now the thermal ratio of Pi to P 2 , and the ratio of P 2 to P 3 , is broken. 

Whilst the equations with three dS vacua are more complicated than the case for two, it 
now appears that the interpretation can be straightforwardly extended. The ratio P 2 /Pi is 
not affected (to leading order) by the presence of the extra minimum, and again its form is 
prescribed by the magnitude of r ls relative to r 2i . The new ratio P 3 /P 2 is found to take its 
thermal value if P2/P1 has a thermal ratio. If P2/P1 has its thermal ratio broken then the rate 
Ti 2 becomes negligible, and the minimum 1 acts as a sink for the remaining two dS vacua. The 
ratio P 3 /P 2 can then be calculated as if they were a system of two dS spaces with vacuum 2 
being allowed to decay to a sink. 

5.3 A Simple "Multidimensional" Landscape 

We now wish to consider the case of a multi-dimensional potential, which we expect to be a 
slightly more realistic model of the landscape. Again, we work with an extension of the simple 
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case of two dS minima and one sink that was given in section [3j Now we consider the potential 
shown in Fig. [51 Here transitions are allowed between the minima labeled 2 and 3 and the 
lower minima 1, which is allowed to decay directly to the sink. This setup models a simple 
potential with more than one dimension, i.e. with more than just pairwise connections between 
vacua. The equations for this potential are given, and solved, in Appendix A. 2. 2. 




Figure 5: A multi- dimensional potential with three positive and one negative minima 
Under the reasonable assumptions that T 2 i Ti 2 and T 3 i ^> Ti 3 , we find: 

P 2 Pi 

— ~ — < 1 and — ~ — < 1 

PL 1 21 "\ 1 31 

when T 2 i and T 3 i ^> Y\ s \ whilst for T\ s and r 3 i T 2 i we find: 

P2 Fig p 3 r 13 

— ~ — > 1 and — ~ — < 1. 

-t"l 1 21 M 1 31 

The case Yi s and T 2 i r 3 i can be found by symmetry from the above expressions, under the 
transcription 2^3. 

These results have a straightforward, but slightly counter-intuitive, interpretation. When 
the sink is narrow with respect to the transition rate T 21 , as well as the rate r 31 , the thermal 
ratios of both P 2 and P 3 , with respect to Pi, are maintained, as expected. When the sink is 
wide with respect to one of the Ts (and narrow with respect to the other), one of P 2 and P 3 
will maintain its thermal ratio whilst the other does not, again, as expected. One could expect 
that when the sink is wide with respect to all other transition rates, all dS vacua will be out of 
equilibrium. Surprisingly, we have found that even in this case only one of the upper vacua is 
out of its thermal ratio with P 1; and the other is not. Furthermore, the vacuum that is taken 
out of thermal equilibrium with P x is the one with the slowest rate of decay. 

We interpret this in the following way: If the sink is wide with respect to all other transition 
rates, then the magnitude of P in the vacuum that decays the quickest will quickly become very 
small. The slowest decaying vacuum, despite the fact that the rate T is smaller, will then be the 



14 



source of the majority of the probability flux at late-times. The small flux that is required to 
keep the faster decaying vacuum in its thermal ratio with the lower vacuum is then maintained 
by the flux from the slower decaying vacuum, which is forced out of its thermal ratio by the 
sink. 

The picture of the slowest decaying vacuum dominating the behavior of the universe may 
seem counter intuitive, but in fact it has a very simple interpretation: Those who want to 
survive in a desert (i.e. near wide sinks) should save water. 

Numerical simulations show that this behavior extends to i > 2 vacua decaying into a single 
lower vacuum, which decays to a sink. If the sink is narrow with respect to all of the other 
transition rates, then the asymptotic ratios P%/P\ approach the thermal ratios that they would 
have in the absence of the sink. If the sink is wide with respect to any or all of the transition 
rates to the higher vacua, then it is only the slowest decaying vacuum that is forced out of 
its thermal equilibrium with the lower vacuum. In this case, the out-of-equilibrium slowest 
decaying vacuum will feed all other vacua, which will be in a state of thermal equilibrium with 
each other. 



We now have some ideas on how to extend the simple model of two dS spaces and one AdS 
space to more general models with multiple dimensions and higher vacua. However, so far 
we have only studied potentials in which one vacuum is allowed to decay to AdS space. In a 
realistic model of the landscape it is likely that many vacua will be allowed to decay in this 
way, and it is the effects of this which we now study. 

5.4.1 2 dS and 2 AdS 

The potential with two dS minima, shown in Fig. [2J can be simply extended to allow both of 
the dS minima to decay to sinks. This situation is solved in Appendix A. 2. 3. Again, we assume 

ri2 < r 21 . 

It can now be shown that if both the sinks are narrow (Tj S <C T^) then the ratio P2/P1 
maintains its thermal value, to leading order. Similarly, if one sink is narrow, and the other is 
wide (Tj S ^> Tij), then the leading order terms are the same as in the absence of the narrow 
sink. Now, if both sinks are wide then the ratio P2/P1 is determined by the wider sink. For 
example, if we have the hierarchy T ls ^> r 2s 3> r 2 i, then 



5.4 Multiple Sinks 



P2 
Pi 



r 2 



21 



.s 



whilst for T 2s ^> T ls ^> T 21 , we have 



P2 r 
Pi ~ r 
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5.4.2 3 dS and 2 AdS 



We may also consider the case of a single higher vacuum that can decay to multiple lower 
vacua, which are, in turn, able to subsequently decay to sinks. We will model this situation 
by considering a potential such as that shown in Fig. [61 where the vacua labeled 1 and 3 are 
allowed to decay to sinks. The equations for this potential are investigated in Appendix A. 2. 4. 
We now have the heirarchy T21 3> V\2 and r 23 3> r 32 . 




Figure 6: A potential with three positive minima, the lower two of which can decay to sinks. 
When both sinks are narrow, Ti s <C T 2 i and T 3s <C r 23 , the thermal ratios 
4 = # » 1 and ^ = f- » 1 

r 2 1 12 * 2 1 32 

are maintained. For one narrow sink and one wide (r^ 3> r 2i ) the narrow sink becomes 
irrelevant, and the problem reduces to the extended one-dimensional potential considered above. 

For two wide sinks, if we take r ls 3> T 3s , without loss of generality, then we always have 
the ratio 

^ = ^«i. 

Pi r t . 

The second ratio, P3/P2, is then determined by the relative magnitudes of r 3s and r 2 i- If 
^3s r 2 i, then we have 

P2 r 3s * ij 

whilst for r 3s <C T 21 we have 

Pr Tqi 

— = — > 1. 
p% 1^32 

This can be understood in terms of the previous example of two dS spaces, both connected 
to sinks. If the probability charge in vacuum 1 is rapidly depleted by the fast decay rate to 
the sink, T ls , then vacuum 1 will subsequently act as a sink for the remaining two vacua. If 
decays to this new effective sink, T 21 , are faster than the decay rate T 3s , then this new sink will 
dominate. Otherwise, T 3s will dominate. 
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5.4.3 Islands 



Consider the potential shown in Fig. [7J Here there are two dS minima which are allowed to 
decay to a sink. (Each of the lower dS vacua could be considered to decay to different sinks, or 
to a common sink, the results are the same). We have also included two higher vacua connected 
to the lower ones, and transitions are allowed between the two lower minima. Here we hope to 
see the effects of having more than one vacuum decaying to a sink. We expect that the presence 
of the second sink will complicate matters and allow for the possibility of two 'islands' that are 
in thermal equilibrium internally, but not with each other. 




s 

Figure 7: A potential with two positive minima that are allowed to decay to anti-de Sitter space. 

With our knowledge of the general solution, we know that the asymptotic attractor solutions 
will be of the form: 

P2 = c ie - mt } 

P 3 = nCl e- mt = nP 2 . 

As shown below, under the reasonable assumption that r 3 4 <C and r 2 i <C Ti2, we find 
three possible values of n and m in Appendix A. 2. 5. In order to choose which solution for 
n = P3/P2 is the relevant one, for a given set of Ts, we must evaluate which m is the smallest. 
This m will correspond to the slowest decaying mode, which dominates in the limit t — > 00; 
the corresponding n will then give the appropriate asymptotic ratio of P 2 /P 3 . A flow chart is 
given in Appendix [Bj which can be used to quickly identify the relevant asymptotic form of n 
for any given set of Ts. 

From the flow chart it can be seen that a necessary condition to maintain a thermal ratio 
between vacua 2 and 3 is that either r 32 ^> T 3s or T 23 ^> T 2s . A sufficient condition for 
maintaining this ratio is either r 23 or T 32 ^> T 2s and T 3s . Conversely, a necessary condition 
for breaking the thermal ratio between 2 and 3 is either T 2s or T 3s ^> T 23 and T 32 ; whilst a 
sufficient condition is given by T 2s and T 3s ^> T 23 and r 32 . 
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Two sets of solutions are relatively simple: 

T 2 3 



ni = 



n 2 = 



F32 + r 3s 

1^23 + ^2s 



F32 

with the corresponding values of m being given by 

mi = r 12 , 
m 2 = r 43 . 

The expressions for m 3 and n 3 are somewhat more complicated, but they can be simplified 
depending on the relative magnitudes of r 23 , T 32 , T 2s and T 3s . If T 23 or T 32 are the fastest of 
these rates then the leading order contributions to m 3 and n 3 are 

T23 , ^2s^32 + ^1^23 

n 3a — — and m 3a ~ 



if T 2s is the fastest then 



T32 

and, similarly, if T 3s is the fastest then 

T23 



r32 r23 + 1^32 



n 3b ~ ^ and m 3b ~ T 32 + T 3s 



c - -p and m 3c ~ T 23 + T 2s ■ 

J- 3s 



We will now consider the significance of each of these three solutions. We take solutions 
rii and mi as corresponding to the situation in which the probability flux is sourced by the 
slowly decaying vacuum 1; solutions n 2 and m 2 correspond to the situation in which vacuum 
4 is the slowest decaying, and therefore sources the probability flux; and solutions n 3 and m 3 
correspond to the situation where the majority of the probability flux is from the lower vacua, 

2 and 3. This interpretation is supported by the various values of m. For m 1 and m 2 it can be 
seen immediately that the asymptotic rate of decay of the Pj is prescribed by the rate of decay 
from 1 to 2 or from 4 to 3, respectively. This strongly suggests that the vacuum dynamics of 
the systems corresponding to these solutions are dominated by the flux out of vacua 1 and 4. 

The interpretation of the m 3 /n 3 solution is a little less straightforward due to its more 
complicated form. When the rate of decay T 2s is fast, m 3 b is the relevant solution. It can be 
seen that m 3 b is determined by the probability flux out of vacuum 3 (recall that r 34 has been 
neglected). By symmetry, when r 3s is large m 3c is the relevant solution, which is given by the 
rate of decay out of vacuum 2. When both rates of decay to the sink are small, m 3a is the 
relevant solution, which can be seen to be some weighted flux out of both of the vacua 2 and 

3 (analogous to the 'reduced mass' of two body dynamics). This justifies the above statement 
that the solutions m 3 and n 3 correspond to a system which is dominated by the flux out of 
vacua 2 and/or 3. 
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We again see that the asymptotic form of the 'vacuum dynamics' is not prescribed by the 
fastest decaying vacua, as may have been naively expected, but by the slowest. These vacua 
hold the majority of the comoving volume at late times, and so dominate the late time evolution 
of the universe. 

The late-time evolution of P 2 and P3 straightforwardly gives the behavior of P± and P4. 
Using the notation p = P1/P2 and r = P4/P3, we obtain 

/ 1 r2i 

Pi t Pi = =r~ and P3 = 

1 12 1 12 1 12 

and 

T34 , 1^34 j 1^34 

ri = ~^ j r 2 =£ — and r 3 = — , 
J- 43 1 43 1 43 

where subscripts i denote that the solution corresponds to m, and rij. It can be seen that if 

the probability charge is being held in vacua 1 or 4, then that vacuum is out of its thermal 

ratio with the lower vacuum to which it can decay, and that the two vacua on the other side of 

the sink are in their thermal ratio. If the probability charge is being held in either or both of 

the lower vacua, then both of the higher vacua have thermal ratios with their respective lower 

vacua. 

Since we have here a large variety of possibilities, let us single out some of the most inter- 
esting regimes. Suppose first that r 23 = r 32 = 0, as discussed in Section HI In this case, in the 
wide sink regime, vacua 1 and 2 will be out of thermal equilibrium with each other, vacua 3 
and 4 will be out of thermal equilibrium with each other, and branches (1,2) and (3,4) will be 
totally disconnected. 

Now let us establish some contact between these branches. If r 23 and r 32 are sufficiently 
small, we will have branches (1,2) and (3,4) out of thermal equilibrium with each other. If 
vacuum 1 has the slowest decay rate among all vacua, it will be out of equilibrium with vacuum 
2. However, in this case vacuum 4 will be in thermal equilibrium with vacuum 3, even if T 2 3 and 
T3 2 are extremely small. Note that the transition to vanishing T 2 3 and is discontinuous. This 
paradoxical situation is similar to the one encountered earlier: The slowest decaying vacuum 
dominates the evolution of other domains, but this vacuum becomes irrelevant if it totally 
decouples from other vacua. Other possibilities can be read from the flow chart in Appendix 

El 

The model discussed in this section shows explicitly that tunneling to sinks can break the 
detailed balance between vacua that would otherwise keep different regions of the landscape 
in thermal contact. This provides a mechanism by which the landscape can be separated into 
different thermally isolated parts, each with a different 'temperature'. In a landscape with 
high enough dimensionality there may be more than one decay channel between 2 and 3. This 
will make it more likely that thermal contact will be maintained between various parts of the 
landscape. 
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6 Discussion 



We have, in the course of this work, carried out a detailed exploration of some toy models of 
the landscape, and studied their thermal properties. Whilst these models are somewhat limited 
in their scope, they have revealed to us a rich, and sometimes counter-intuitive structure. 

In general, it appears that a necessary condition for the disruption of thermal ratios between 
vacua is a rapid rate of decay to an AdS space - a wide sink. Intriguingly, however, once 
this condition is met, the late-time behavior of the system is controlled by the slowest of 
the tunneling processes. For example, when we considered a simplified version of a multi- 
dimensional landscape with a single sink (Section 15. 3p we found that only the slowest decaying 
vacuum was shifted from its thermal ratio with the lower vacuum: All other ratios remained in 
thermal equilibrium. We interpret this as being due to the more rapidly decaying vacua quickly 
depleting their populations, leaving the slowest decaying vacuum to be the primary source of 
probability flux. The late-time thermal behavior of the multiverse is then determined by the 
weakest transitions - transitions that would otherwise be irrelevant in the absence of a sink, or 
if we took the limit of their rate going to zero. 

We see a similarly interesting interplay in our "islands" example in Section l5.4.3l Here, once 
thermal equilibrium is disrupted by a wide sink, the late-time behavior of the system is again 
determined by the slowest decaying vacuum. In this example the sink can break the thermal 
contact between different regions, leaving several thermally disconnected 'islands'. 

In a more complex model there may exist intermediate vacua that can act as bridges, 
restoring thermal contact between islands. Such vacua would be required to be able to decay 
to both of the islands, with a non-negligible rate. As we have shown, quickly decaying vacua 
are likely to be in thermal equilibrium with the lower vacua to which they are allowed to decay. 
It may happen that for a realistic landscape with many dimensions and very many vacua, the 
thermally isolated regions will be relatively rare. However, in order to verify this conjecture 
one would need to study simultaneously a disruptive effect of a very large number of sinks and 
a restoring effect of a very large number of bridges. 

We anticipate that further exploration along the lines of this paper may reveal still more 
interesting features, allowing us a more complete and better understood picture of the land- 
scape. 
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A Solving the Equations 

In this appendix we will elaborate on some mathematical details of the results given in the 
main body of the text. 
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A.l General Solution 



As discussed in section El the equations governing the evolution of some set of Ps are: 

§ = -Er«« + Eivi-™. 

3Ti 3+i 

The first term on the right hand side gives the flux out of the state denoted by Pi to all the 
other Pj, the second gives the flux into Pi and the third gives the flux out of Pi and into the 
sink. For n different Ps, we have a set of n coupled first order ordinary differential equations. 
Such a set of equations can be manipulated into a single linear nth order equation, for any one 
particular Pj, of the form: 

£%^r = °- (20) 

3 

The coefficients aj can be written in terms of the transition rates, T. This equation has the 
general solution 

n 

I', • (21) 

3 

where the Cj are constants of integration and the constants rrij are the n roots of the nth order 
polynomial YTj djm J = 0. All other Pj can then be seen to have the same functional form, by 
substitution back into the original set of equations (fT8j) . The coefficients of the different modes 
of these other Pj will be completely determined in terms of Cj and the Ts. This gives the general 
solution to (118I) . for all Pj, with n arbitrary constants. 

At late times the dominant mode in the solution (|21|) will be the one with the smallest 
exponent, m. As all the Ps have the same functional forms, the same mode will dominate the 
evolution of each P at late times, and so the ratios Pi/Pj will asymptotically approach constant 
values. The form of the dominant mode can now be calculated in terms of the constant ratios 
qij = Pi/Pj by adding all the equations (|T8|) to get 



:22) 



i i 

which can be integrated to 

InP, = - Sk±S^) (t _ £o) , (23) 

where to is an arbitrary constant. We now see that the smallest value of m, which corresponds 
to the dominant term as t — > oo, is given by 

m = ^; + ^ r ^) < o. (24) 
Therefore the real part of all of the roots m, are > 0, and all the modes are decaying. 
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We now know that the late-time attractor solution for each of the Ps is an exponentially 
decaying function, and that this rate of decay is the same for each P. This allows us to work 
out the constant ratios g^, that are asymptotically approached as t — > oo. Taking the mode 
with the smallest m and substituting it back into (|T8|) gives us a set of n — 1 algebraic equations 
which can be solved for the n — 1 unknown ratios g^. The asymptotic form of the Ps is now, 
up to a normalization, completely solved for. 



A. 2 Details of Particular Solutions 

Here we will provide the equations that were used to find the results in the main body of the 
text, above. We will assume that all transition rates, T, are orders of magnitude different from 
each other. 



A. 2.1 An Extended ID Potential 

The equations governing the system shown in Fig. H] are given by 

Pi = -(ri2 + r ls )p 1 + r 21 p 2 

p% — — (T21 + r 23 )p 2 + r 12 Pi + r 32 -P3 

P3 = — F 32 P 3 + F 23 P 2 • 

We know that the solutions will be of exponential form, so we substitute the ansatz 

P 2 = Cl e~ mt 
Pi = pcie~ mt = pP 2 
P 3 =q Cl e- mt =qP 2 , 

and find the expressions 

= r 23 + g(r 23 + r 21 - r 32 - pv 12 ) - v 32 q 2 (25) 

= r 2 i + p(r 2 i - r i2 - r ls + r 23 - gr 32 ) - r 12 p 2 , 

or 

= (r 32 g-r 23 )(l + g+p)-r ls pg (26) 

= (T 12 p-T 21 )(l + q + p)+T ls p(l + q) . 

We find three sets of solutions. Firstly, when T ls either T 12 or T 21 , we see from equations 
(ESD that: 

r 2 i r 23 

p — — — 3> 1 and q ~ — — , 

Ti 2 r 32 

where we have taken T 21 ^> T 12 . In this limit the thermodynamic ratios (Q are maintained. 
For p 1 and g <C 1, p and q are given by (1261) as: 

r 2 i r 21 

p ~ ~ — i 

Ti 2 + r ls Tig 

(Ti 2 + ri s )r 23 r 23 r 23 

q ~ ~ ~ — 1 

Ti 2 r 32 + ri s (r 32 — r 2 i) r 32 — r 2 i r 32 
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where (r± a ^> r 12 and T 2 i) and (r 32 T 2 i and r 23 ). Here the thermal ratio of Pi to P 2 is 
broken, whilst that of P 2 to P3 is maintained. Lastly, for p C 1 and q ^> 1, the solutions to 
(IH are: 

r 21 

p ~ ~ — <c 1 

Ti2 + Tis Tis 

1^121^23 + ri s (r 2 i + r 23 ) r2i + r23 ^ -, 

q — ~ ^> 1 

(r 12 + r ls )r 32 r 32 

where (r ls 3> T 12 and T 21 ) and (r 32 <C either T 21 or T 23 ). Again, the thermal ratio of Pi to P 2 
is broken and now the ratio of P 2 to P 3 is thermal if r 2 i <C T 23 and not thermal if T 2 i ^> r 23 . 
Note, that no solutions were found for p and q ^> 1, when the sink is wide. These results have 
been verified numerically. 



A. 2. 2 A Simple Multi-Dimensional Potential 

For the setup show in Fig. [5] the system of equations governing the vacuum dynamics is 

A = -(ri 2 + ri 3 + ri a )Pi + r 2 iP 2 + r 3 iP 3 



P2 — - T2iP2 + r^Pi 
P3 — - r 3 iP 3 + r^Pi 



Substituting the ansatz 



P x = Cl e- mt 

P 2 = pc\e~ mt = pP] 



gives the equations 



P 



Q 



P3 — QCie = qP\ , 

(r 3 ig-r 13 )(i + g)-gr lfl 

Ti 3 - r 3 ig 

(r 2lP -r 12 )(i+p)- P r ls 



r 12 - t 2 ip 

Under the reasonable assumptions that T 2 i ^> r 12 and T 31 ^> r 13 , for p <C 1 and g < 1 we 
have 

P - p w~ - tt~ < 1 and q ~ — ~ — < 1 

1 21 ~ 1 Is 1 21 1 31 ~~ 1 Is 1 31 

when T 21 and T 31 3> T ls ; whilst for p 3> 1 ^> q we find: 

Ti2 + Tis Tis , 

P - — p p— > 1 

J- 21 1 21 

ri 3 (ri 2 + r la ) r 13 

g ~ ~ <C 1 

1^121^31 + risi^si ~~ ii s r 2 i r 3 i 

when Tis and T 3 i ^> r 2 i. The case q ^> 1 ^> p can be found by symmetry from the above 
expressions, under the transcription 2<-^3 and p <-> q. We find no consistent solution exists with 
p> 1 and g > 1. These solutions have been numerically verified in the limits indicated above. 
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A. 2. 3 Two de Sitter and Two Sinks 



The dynamical equations for this situation are given by: 

Pi = -(r 12 + r ls )Px + r 21 p 2 
p 2 = -(r 2 i + r 2s )p 2 + r 12 P! , 

which have the general solution: 

P2 r 12 - T 21 + T ls -T 2a + B tanh [f (t - t )] 

Pi 2T 21 

_ r i2 — r 2i + Vis — r 2s + b 
~* 2r 21 

as t — > 00. Here B 2 = (r i2 + T 2 i + T u + r 2s ) 2 - 4(Fi 2 r 2s + Ti s r 2 i + Fi s r 2s )) and the results 
in the text follow from this asymptotic form. 

A. 2. 4 Three de Sitter and Two Sinks 

The relevant evolution equations are now: 

A = -(ri 2 + r ls )P! + r 21 p 2 

Pi — — (T21 + r 23 )p 2 + r^Pi + r 32 p 3 

P3 = ~(F 32 + F 3s )P 3 + T23P2 



which yield 



= r23 + g(r23 + r 21 -r 3s -r32-pr 1 2)-r32g 2 (27) 
= r 2i + p(r 2 i - r 12 - r ls + r 23 - <?r 32 ) - r 12 p 2 , 



or 



(l + q)( P Ti s + qT 3s ) = (1 + p + q) {qT 3s - pT u + T 21 ) (28) 
(l+p)(pTi s + qT 3s ) = (l+p + q )(pTi s -qT 32 + T 23 ) 

where p = Pi/P 2 and q = P 3 /P 2 . When both sinks are narrow (ri s <C Ti2 and T 3s <C r 3 2) it 
can be seen from ( |271) that the thermal ratios 

£ = and J = ^»l 

r 2 1 12 * 2 1 32 

are maintained. For one narrow sink and one wide (T is <C T i2 ) the narrow sink becomes 
irrelevant, and the problem reduces to the one considered above. 

The case of two wide sinks now remains. From (|28p . for p <C 1 and g C 1 we obtain: 

T2i T 2 i 
p ~ ~ <C 1 

Ti2 + Tis r ls 

T23 T 2 3 

g ~ ~ <C 1 

r 3 2 + r 3s r 3s 
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whilst p <C 1 <C q gives: 



r 2 i T 2 i 
p ~ ~ <C 1 

Ti2 + ii s Ti s 

T23 + ^21 — ^3s T 2 l 

g ~ ~ >1 

-L 32 J- 32 

where r 32 <C r 23 <C T 3s <C T 2 i <C T ls . The remaining case p ^> q ^> 1 has no solutions (given 
two wide sinks). 



A. 2. 5 Islands 

The potential shown in Fig. [7] is governed by the set of equations: 

-Pi = — Tl2-Pl + T21-F2 

Pi — — (T21 + r 23 + r 2s )-P 2 + r 12 Pi + r 32 p 3 

P3 — — (T34 + + TssjPs + T^P^ + r 2 3-P2 

Pa — — F43P4 + T 34 P 3 . 

These four coupled first-order equations can be recast into a set of two coupled second-order 
ordinary differential equations, for the variables P 2 and P 3 

P2 + (T12 + r 21 + r 23 + r 2s )p 2 + r 12 (r 23 + r 2s )p 2 = r 32 p 3 + r 12 r 32 p 3 
P3 + (r43 + r34 + r32 + r 3 ,;)p 3 + r 43 (r 3 2 4- r 3s )p 3 = r 23 p 2 + T4 3 r 23 p 2 . 

We know that the asymptotic attractor solutions have the form 

P 2 = c ie - mt 

P 3 = n Cl e- mt = nP 2 , 

which on substitution into the second-order equations above, give 

r i2 (r 23 + r 2s - ™r 32 - m) - m(r 2 i + r 23 + r 2s - nr 32 - m) = o 
n((r 32 + r 3s )(r 43 -m) -m(r 3 4 + r 43 -m))-r 23 (r 43 -m) = o. 

Under the reasonable assumption that T 34 <C T 43 and T 21 <C T 12 these two equations can be 
solved to give the three possible values of n and m discussed above: 

r23 

n x -- 



n 2 



1^32 + T3s — Tl2 
1^23 + ^2s ~ ^43 



^32 

^23 — ^32 + ^2s — ^3s + a/ (^23 ~ T 32 + T 2s ~ T 3s ) 2 + 4r 23 F 3 2 



2T 32 
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with the corresponding values of m being given by 



mi = Ti2 

m 2 = r 43 

^3 = 2 ( r 23 + r 32 + r 2s + t 3s - a/ (r 23 - r 32 + r 2s - r 3s ) 2 + 4r 23 r 32 ) . 

There is a fourth mathematically permissible value of m and n, however it corresponds torn < 
for all Tij and so we do not consider it to be of any physical significance (a negative value of P 
makes very little sense). These results have been confirmed numerically. 
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B Flow Chart 



Fig. [8] shows the various possible asymptotic limits for the "thermal islands" discussed in 
section 15.4.31 



Which of 

(a) r 2 3/r 3 2, (b) r 25 or (c) r 3S 

is largest? 




or c 



Which of 

» r 12 , (b)r 43 

! r 2 gr32+r 3 5r23 
is smallest? 



(b) 






Which of 

(a) r 12 , (b)r 43 
or (c) r 25 + r 23 

is smallest? 



Which of 

(a) r 12 , (b)r 43 
or (c) r 35 + r 32 

is smallest? 




T23 

Ul ~ r 3S 



m „ r ~^~ / -^23 

712 ~ ril 






n-i ps 
° r 3S 



n 2 w r 



r2 3 +r 2 g 



n 2 



r 2 g 

T32 



Figure 8: Flowchart showing different asymptotic limits for the "islands" example. 
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